QUASILIMITING BEHAVIOR FOR ONE-DIMENSIONAL DIFFUSIONS 

WITH KILLING 



MARTIN KOLB AND DAVID STEINSALTZ 

Abstract. This paper extends and clarifies results of Steinsaltz and Evans [47], which found 
conditions for convergence of a killed one-dimensional diffusion conditioned on survival, to a 
quasistationary distribution whose density is given by the top eigenfunction of the generator. 
Convergence occurs when the limit of the killing at cxd differs from the negative top of 
the spectrum of the generator. When the killing at oo is larger than the negative top of 
the spectrum, then the eigenfunction is integrable. When the killing at oo is smaller, the 
eigenfunction is integrable only when the unkilled process is recurrent; otherwise, the process 
conditioned on survival converges to density on any bounded interval. 



1. Introduction 

1.1. Background and history. Killed Markov processes are central objects in probability 
theory. When considering their asymptotic behavior there are two distinct approaches, which 
may be thought of as two different ways of ordering the different limits: 

• We condition on the process never being killed. That is, we look at the distribution 
of {Xf : t £ [0, s]} for fixed s conditioned on tq > T (where tq is the killing time), in 
the limit as T — )■ oo; we may then take a second limit s — )■ oo to define the process on 
[0, oo). Following Doob [IT], this is called an h-process or an h-transform of the 
original Markov process, and it generally produces a new Markov process. A well- 
known example of this is the 3-dimensional Bessel process, which may be derived 
from the one-dimensional Brownian motion, conditioned never to hit [50' Section 
6.6]. This is intimately connected to questions about the Martin boundary. 

• We condition on the process not having been killed up to the time in question. That 
is, we look at the distribution of Xt conditioned on rg > t, in the limit as t — )• oo. 
The collection of distributions for different times t are not consistent, and so cannot 
be analyzed directly with Markov-process techniques, but are more amenable to an 
analytic semigroup approach. We generally expect the distribution to converge to 
a limit — sometimes called the Yaglom limit, after the seminal work of Yaglom 
|56] on branching Markov processes conditioned on long survival — which will be 
quasistationary, in the sense that when started in this distribution the process will 
remain in a multiple of the same distribution for all times. 

The extensive mathematical development and wide-ranging applications in this 
bibliography of papers on quasistationary distributions and Yaglom limits compiled and pe- 
riodically updated by Phil Pollett |42) lists 375 entries through 2005 — permit us to mention 
only a smattering of the vast array of applications of killed Markov processes to biology 



2000 Mathematics Subject Classification. Primary 58J05; Secondary 35J10, 35P05, 47F05. 
Key words and phrases, diffusion, killing, quasistationary distribution, spectral theory, Sturm-Liouville, 
operator. 



2 



MARTIN KOLB AND DAVID STEINSALTZ 



[g [231 [221 [9], demography [Ml [291 [30], medicine [Ml [5Z|, and statistics P [3 [l5]. Par- 
ticularly in the demographic and medical contexts, where killed Markov processes suggest 
themselves as models for populations undergoing culling by mortality or other processes, Ya- 
glom limits, while rarely mentioned explicitly in the applied literature, correspond naturally 
to the observable distribution of survivors. 

The central concerns of this theory are to describe, for a given class of sub-Markov processes, 
the quasistationary distributions (if any), and to describe the convergence (or not) of the 
process conditioned on survival to one of these quasistationary distributions. A significant 
part of the literature focuses on discrete state spaces, commonly birth-death processes (or 
with some more flexible localization of the transitions), with killing only on the boundary. 
|18] gave one of the most general accounts of the existence and convergence to quasistationary 
distributions for discrete processes of this kind. One unusual contribution, outside of these 
categories, is [21], which proves convergence to quasistationarity for fairly arbitrary discrete 
Markov chains with general killing, by imposing a stringent Lyapunov-like drift condition. 
|28] generalize in a different direction, finding existence and vague-convergence conditions for 
discrete-time Markov chains on general metric spaces. 

Killed birth-death processes naturally generalise in the continuous-space context to killed 
diffusions, which have received rather less attention. The existence of eigenfunctions for the 
generator of a one-dimensional diffusions is simplified in the continuous setting, as we may rely 
upon standard theory of ordinary differential equations. Showing that these eigenfunctions are 
integrable (hence represent the densities of distributions) and quasistationary is more involved, 
though, and showing Yaglom convergence to the minimal quasistationary distribution becomes 
technically challenging, particularly when the state space is an unbounded interval. The 
basis for all later work on Yaglom convergence of diffusions was laid by Mandl [33], who 
used standard results from Sturm-Liouville theory and the spectral theorem for self-adjoint 
operators to prove vague convergence (that is, convergence of the distribution of the process 
conditioned on being in a compact set), and uniform convergence under an assumption of 
strong inward drift. These results have been substantially extended by a shifting coalition of 
researchers who have produced papers 112 ^ 135 ^ 19]. which elucidate the conditions under which 
Yaglom convergence occurs, and distinguish in [36| the R-positive situation from the R-null 
— essentially, exponential-rate decline of probabilities distinguished from declines that are 
asymptotically not exactly exponential — in terms of the eigenfunctions. 

One important constraint in nearly all of this work — as well as Pinsky's results in [39], 
for diffusions on a compact domain with gradient-type drift — is the assumption that killing 
occurs only at the boundary. Not only is this an unnatural assumption from the perspective of 
many of the applications, particularly the demographic applications, as discussed in [iTj [36]. 
we argue here that it obscures the fundamental links among the spectrum, the killing rate 
out at infinity, the recurrence-transience dichotomy, and Yaglom convergence. (An exception 
which proves the rule is the biological application of internally killed diffusions |25J, which 
makes no reference to any of the literature on killed diffusions, and cites only [l3] for qua- 
sistationary distributions of discrete chains.) In |47j it was shown that a key criterion for 
Yaglom convergence is whether the limit of the killing rate out at infinity is greater than or 
less than the bottom of the diffusion operator spectrum Aq, which may be interpreted as the 
rate of being killed or escaping from a compact interval around the origin. Here we substan- 
tially complete that program of |47J, providing testable criteria for distinguishing convergent 
processes from divergent processes that escape to infinity when conditioned on long survival, 
in all cases when the limit killing rate at infinity exists and is not exactly equal to Aq. 
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1.2. Heuristics. We begin by observing that general spectral theory — summarized here in 
Lemma [2?T] — tells us that the bottom of the spectrum Ag gives the exponential rate of decay 
of the distribution of Xf restricted to a compact interval. What needs to be addressed, then, is 
the question of whether the portion of the surviving mass within a compact interval dominates 
the total surviving mass. There are two ways of addressing this question. One is in terms 
of the spectrum of the 2'^ generator. Consider the case when Ag is an isolated eigenvalue, as 
when the diffusion takes place on a compact interval with two regular boundaries, so that the 
resolvent is compact; but also in cases which intuitively seem well-approximated by a compact 
process, as when there is strong drift pulling the process in from oo, and when there is strong 
killing out toward oo. Then we expect the same methods that work in finite dimensions, 
for powers of positive symmetric matrices, to work here as well. One catch is that the £p 
convergence need not tell us about the convergence of the conditioned density, which is an 
2^ property; indeed, it is easy to see (c/. Proposition 2.3 of [17]) that the process always 
escapes to oo if ¥''^0 > x)dT{x) is not finite, where T is the speed measure and if{X, •) is the 
eigenfunction of the generator, defined as the solution to an ordinary differential equation in 
section 13. 2i We show, in Lemma 14.21 that these conditions suffice: That is, whenever Ag is 
an isolated eigenvalue, and the corresponding eigenfunction is also integrable, then we have 
convergence to the quasistationary distribution given by the density f{XQ, •)/ f fi^Q, x)dr{x). 

What about the case when Ag is not an isolated eigenvalue? This corresponds to the i?-null 
and i2-transient cases in Tweedie's theory [491 148j . where the decline of the transition kernel is 
not exactly exponential with rate — Ag , but slightly faster, in the sense that e^o*p(i, x, y) — )• 0. 
When Ag is greater than the limit of killing at 00 it turns out that the convergence lines 
up precisely with the standard recurrence/transience dichotomy for the unkilled process. 
(Another way of putting this is to say that when the i2-recurrence or i?-transience does not 
conform to the properties of the unkilled process, this must be reflected in the equality of 
Ag and \\m.x^ao i^[x) .) As summarized in Lemma [3.31 if Ag > lim sup2._i.00 — if the 
distribution declines on compact sets at a faster exponential rate than would keep pace with 
the killing out toward 00 — this implies an exponential decline on compact sets even in the 
absence of internal killing (by which we mean killing by k, rather than instantaneous killing 
at the boundary). This means either that the scale is finite, so that the unkilled process 
runs off toward 00, and the distribution on compact sets would decline exponentially in the 
absence of any killing at all — the i?-transient case — or the scale is infinite with finite speed, 
which means that the process keeps returning to 0, at long intervals, and the killing rate 
Ao corresponds to real killing at 0, not escape, which is the -R-null case. In the i?-transient 
case the conditioned process escapes to infinity. In the i?-null case the conditioned process 
converges to the quasistationary distribution. 

1.3. Main results. The core of this work is the identification of the asymptotic behavior of 
killed diffusions in terms of the eigenvalues of the generator and the finiteness of the speed 
and scale measures. We move beyond earlier work in eliminating the supplementary bounds 
on the drift and killing terms, and in providing easily testable criteria for determining whether 
the conditioned process converges to a quasistationary distribution for all cases in which the 
bottom of the spectrum Ag does not coincide with the limit of the killing rate at 00. 

We begin by summarising the most important results. These results presuppose general 
assumptions and restrictions on the processes involved, which will be formulated fully in 
section [2j The quasistationary distribution will be defined in terms of its density (/?(Ag, •) 
with respect to F, where ip is the eigenfunction of the generator, defined as the solution to an 
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ordinary differential equation with appropriate boundary condition, stated formally in section 
[321 

(i) There is always convergence to the quasistationary distribution on compact sets, 
stated formally as Theorem 13. 151 

(ii) If Aq > lim sup2._^Qj3 k{x) or Ag < lim inf^._!,oo ^(2;), then the conditioned process either 
converges to the quasistationary distribution with density v'(Aq , •) [J^ ip{XQ,y)dT{y)) ^ 
with respect to T, or escapes to 00. This behavior is independent of the initial dis- 
tribution, provided only that it is compactly supported. (For explanation of the 
terminology, see section [2^31 ) This is Theorem 3.3 of [17], but it is restated here as 
Theorem 12.61 in a slightly stronger form, as several restrictions have been removed. 

(iii) If Aq < limsup^^oo ^(x) then the conditioned process converges to the quasistation- 
ary distribution. This is stated as Theorem 14. 1[ Note that this includes the (some- 
what unintuitive) fact that a bound on the £^ spectrum implies that ux^ {z)^{z)dz < 
00, which is a fact about the £^ spectrum. 

(iv) If K := lima;_j.oo k(x) exists and K < Xq, then the behavior of the conditioned process 
depends on the transience or recurrence of the unkilled process. If the unkilled process 
is transient — that is, if ^{x)^^dx < 00 — then the conditioned process escapes 
to 00. If the unkilled process is recurrent — that is, if ^{x)~^dx = 00 — then 
the conditioned process converges to the quasistationary distribution. These results 
are stated in Theorems 14.91 and 14.71 

These extend similar results from [17]. The differences are as follows: 

• The boundary at 00 was assumed in the earlier work to be natural. Here, it may be 
an entrance boundary as well. 

• In [17] there were conditions imposed that required the drift and killing not to grow 
too quickly, or be too irregular. Here there is no constraint on the killing other than 
local boundedness, and no constraint on the drift other than that which implies that 
00 is inaccessible. 

• The earlier work imposed an unintuitive condition on the initial distribution that 
could be hard to check. Here the results all presuppose only a compactly supported 
initial distribution. 

• The most substantial advance: The case K < Xq is now shown to be split by the 
standard recurrence-transience dichotomy, which tells us whether the conditioned 
process converges or escapes. In [47] this was left open. 

2. Assumptions, Definitions and Previous Results 

2.1. Analytic terminology. In general a Sturm-Liouville operator is any formal differential 
operator of the form r = Tp^qy = ~^'h'ilx + ^' "^^ere p, q : {a, b) —?■ (0, 00) and V : {a, b) — > 
M are sufficiently well behaved functions. In this work we consider only operators where 
p = q = 'j,V = K>0 and a = 0, 6 = 00. Moreover, we always assume in this chapter 
that 7(x) = e^/o'^'W'^^ for some b G 2l^{[0,oo)) D C((0,oo)) and < k E C([0,oo)). These 
conditions are not entirely necessary, but this constraint still admits a large class of one- 
dimensional diffusions. (However, see [9] for a natural application to biology which requires 
b to be singular at 0.) Concerning the assumptions on b we could replace the condition 
b G £,l^^{[0,oo)) by the condition that Jq e~ ■fc '^^^^^''■^ dx < 00 for some c G (0, 00), and 
< 00, which is equivalent to saying that the boundary point is regular in 
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the sense of Feller and also in the sense of Weyl. The formal differential operator L'^ 
— ^^7^ + K gives rise to a closable densely defined quadratic form in £^ by 

jam' + \ /o~ W\y)?l{y) dy + /o°° <ymy)? i{y) dy if a < oo, 



h ir l'/^'(y)P7(y) dy + K(y)|99(y)p 7(y) dy if a = oo. 



(1) yp^r-'M 

for any £ V^, where a = po/{l — po) £ [0, oo], and the domain is 

T^a '■= {y? S £^ I (f, jif' absolutely continuous in (0, oo), q'^{ip) < oo, (1— po)0(O) = ^Po'P'i^)} 

The closure of this quadratic form will be denoted by g"'". In the sequel T will denote the 
measure 7(7/) dy. (This is identical to the "speed measure" of the diffusion.) Where there 
is no need to emphasise the particular boundary condition at 0, the parameter a will be 
suppressed. 

In this paper we will consistently use F as a reference measure instead of the Lebesgue 
measure, which is different from the convention adopted in [47|. Recall that the speed mea- 
sure of a one-dimensional diffusion is also the reversing measure, with respect to which the 
generator is symmetric. Unless otherwise indicated, we will always use the bare notation £^ 
to mean £^((0, 00), F) , and for f,g £ 2.'^ we have the inner product 

POD 

(2) {f,g) = / fix)gixhix)dx. 

Jo 

To the quadratic form q'^'°' there corresponds a uniquely defined positive selfadjoint operator 
"When a = this is the so-called Friedrichs extension. It is easy to see that the action 
of the operator L'^'" is given by 

L«."(^(x) = -^^"{x) - b{x)^'{x) + K{x)(p{x). 

The bottom of the spectrum of L'^ will be denoted by Aq. The spectrum of the selfadjoint 
operator L'^ is written E(L'^). Where there is no danger of confusion, the corresponding 
objects with k = will also be denoted by q, L and Aq instead of q^, LP and Ag, respectively 
(or gO'", LO''", and Aq'". Since L'^ and L are selfadjoint operators the spectral theorem implies 
the existence of spectral resolutions (-E'^)Ae[Ag,oo) ™d (-E'A)Ae[Ao,oo)) respectively. For the basic 
facts concerning spectral theory of selfadjoint operators the reader should consult [53]. 

The spectral theorem for self-adjoint operators allows us to define functions f{L'^) of the 
operator. For every Borel-measurable function / : M — )• M the operator f{L'^) is defined via 

(3) -DifiL^) = In G £2 I r |_^(^)|2 dWE-ufiX) < 00 



(4) fmu= [ f{X)dE^X)u, 

(5) Vmuf = [ f{Xfd\\E-u\\\X). 

Observe that for a Borel-measurable function / : [0, 00) — )• M and a > we have Ran{f{L'^)) C 
D((L'^)") if [0,00) 9 A I— )• |A"/(A)| is bounded. This implies in particular that the range of 
g-tL'- -g contained in the domain of all powers of L". Moreover the spectral theorem allows 
us to clarify further the connection between the quadratic form q'^ and the associated non- 
negative operator L'^. Let ^/L^ denote the unique non- negative square root of L'^, which is 
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defined using tiie spectral theorem. Then we have 'D(q'^) = and for every / G T>{L'^) 

we have 

(6) q-{f,g) = {VL^f,^g). 

Using the 'elliptic' Harnack inequahty and Weyl's spectral theorem it is not difficult to see 
that 

Aq = max{A G M | there is a positive solution of (L'" — X)u = 

withn(0) = -^, ^'(0) = -^}. 

1 + a 1 + a 

(This was proved by [33] using slightly different methods.) Equation ([7]) already suggests that 
for < A < Ag solutions of {L'^ — X)u = might have a probabilistic significance. 
In the sequel we usually denote by (p{X, •) the solution of the eigenvalue equation 

1 (T 

(8) {L- - XMX, •) = 0, (^(A, 0) = f'iK 0) 



1 + a 1 + a 

It might be important to note that solutions in ([7]) and ([8]) are solutions in the sense of 
the theory of ordinary differential equations. An important issue is whether the solution 
also belongs to the Hilbert space and thus is an eigenfunction in the sense of spectral 
theory. When we wish to emphasise that certain solution are also eigenfunctions in the sense 
of spectral theory, we denote them by ux- 

Crucial to much of our analysis is the fact that the asymptotic behaviour of the semigroup 
is wholly determined by the spectrum right near the base of the spectral measure, which 
we show in Lemma |2.H and then that the base of the spectral measure for any nonnegative 
function is Aq, which is Lemma |2.2[ For g G Q'^, define A^ to be the infimum of the support 
of the spectral measure of g; that is, 

(9) Xg :=sup{A : = O}, 

and let Ax be the subspace of consisting of functions / such that A/ > A. 
Lemma 2.1. Given g G D(L'''"), we have 
\g{x)\<Ca{x)\\VL^g\\+C'M\ 

(10) / roo xl/2 / .oo xl/2 

= Ca{x)^J^ Xd\\E^g\\'{X)j +C'a[l d\\E^g\\\X) 



where 

(11) Ca{x) 
and 



max{y|,(2/o"7(2/)-idy)'/'} fora>0, 
;i8/o"7(2/)-^dy)'^' fora = 0, 



^-^2) C' ■= 12 > ° io°°^(y)^y = 

I iIo° ^iy)dy) '^^^ ifa = and l{y)dy < oo. 

For any t > l/2Xg, 

(13) sup |e-*^"5(2;)| < {Ca{x)Xg + C'^) 115"--*^ 



e 
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Proof. Suppose a G (0,oo). Since g G V{L'^) is differentiable, we have 



\g{x)\ < \gm + / \g'{y) 



MM 
i{y) 



1^(0)1 + 

7 



< |5(0)| + 



7 



< (^2|ff(0)|2 + 2 
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(Cauchy-Schwarz inequality) 

1/2 



\9'{y)\ i{y)dy 



by ([!]) and ([H]). The spectral theorem (0) allows us to represent \rL^g in terms of the spectral 
resolution, yielding ([TO]) . 

If a = oo then g(0) = 0, so the corresponding term drops out of the bound. 

If a = 0, we have the alternative bound 

\g{x)\<\gm^C[f^\g)f' 
\g{x)\>\gm-C{q''\g)i'\ 
. The second bound gives us 

jiy)dy) (|5(0)f-2C|5(0)|(g'^'0(<7))'^' 



where C = V2 



'-I0,x] 



\9f > 



which implies that 



|5(0)| <2C7^^+||5|| 



i{y)dy 



-1/2 



We combine this with the above calculation to obtain the appropriate version of (jlU|) . 
For any positive t, we have gt := e~^^'^ g G 'D{L'^), so we may apply (jlOp to obtain 

1/2 



\9t{x)\ < 



l{y) dy 







-tL- 



Applying again the spectral theorem ^ — now with f{x) = \fxe *^ 

2 



yields 



poo 

/ Xe~^'^d\\E''gf {X) 
Jo 

roo 

= / Ae-2*^d||S^(7f (A) 
<X„e-^'^^\\gf, 



since Ae attains its maximum at A = -^r. Similarly, 



""dWE^gfiX) < e'^^^^\\g\ 
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Lemma 2.2. For any non-negative measurable function / E £^ with \\f\\ > 0, the spectral 
measure d||£"'/|p(A) corresponding to f includes Ag in its support. 

Proof. Since e~^'" f is everywhere nonnegative (except perhaps at the boundary), and its 
associated spectral measure has the same support as we may assume that if there 

were a counterexample it would not vanish, except perhaps at 0. 

Suppose there is some A* > Aq such that = 0. Then for any /i G £^ with \h\ < /, 



> 



> 



-At 



d\\E^ff{X) 



A. 



oo 



AS 



dWE^hnx). 



Thus, it must be that (i||ii^'^/i|p(A) is supported on [A*,oo) as well. Thus, for all such h we 
have \\E^*h\\ = 0. 

Let fn = f • l[n,oo)- For any A G (Aq,A*), by ([3]) is in the domain of the resolvent 
Rj^ = {L'^ — A)^^. Furthermore, by ([5]), if we choose A^,* large enough so that ||Ea,,/|| > 0, 
then 



\R'J\ 



(A - X)-^d\\E-fr{X) > (A - A„)-2||Ea„/|P > 0. 



Let gn := R-^fn/\\Rxfn\\- Then gn satisfies 

L'^9n{x) = Xgn for x <n. 

(In principle, the equality holds only in the S? sense, but it becomes true for all x since both 
sides are in 2?^, hence in particular continuous.) By the representation 



we see that gn is nonnegative. 

The space of solutions to the ordinary differential equation L'^g = Xg satisfying boundary 
condition (I14p is one-dimensional, so if we renormalise to 



9n 



(1 + a) ^fif„(0) ^gn ifa<oo, 
5n(0)"^5'n ifa = oo. 



we have gn{x) = gn'{x) for x S [0,n] when n < n'. The limit must then be identical with the 
function (p{X,-), and is everywhere nonnegative, contradicting the characterisation of Aq in 

□ 
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2.2. Boundary conditions, recurrence, and transience. Defining the diffusion includes 
a boundary condition at 0, parametrised by po G [0, 1] or a = 2(l/po — 1) G [0, oo]: 

(14) (l_po)0(O) = ipo0'(o). 

The condition a = oo corresponds to instantaneous killing at 0, while a = corresponds to 
reflection with no killing. Intermediate parameters correspond to "slow killing" at 0, so that 
the process is killed when the local time at reaches an exponentially distributed random 
variable. The operator L**'" is associated with the closure of the quadratic form q'^'". That 
is, L is the selfadjoint realization of the differential expression (t^) + k in £^ that has 

boundary condition ()14p at 0. The quadratic form g is a Dirichlet form, and the canonically 
associated Markov process is a solution for the martingale problem associated to the operator 
L with the appropriate killing or reflection at 0. This means there exists a family of measures 
(Pt)tg(o,oo) on the space C([0, oo),M) of real valued continuous functions on [0, oo) such that 
for every / G £^ and every x G (0, oo) (due to the Feller property) 

(e-*V)(x)=E,[/(Xt),To>t], 

where (Xt) is the canonical process on C([0,oo),]R) and Tq is a random time defined with 
respect to the local time at 0. (Again, if a = oo then Tq is the time of first hitting 0; if a = 
then Tq = oo.) In this normalisation, the scale measure has density 7(x)~^ with respect to 
Lebesgue measure. 

It is a trivial consequence of the definition of natural scale that J^^ -/{x)^^ is a martingale, 
and so that (^Xt hits eventually) = 1 for x > if and only if the scale function is infinite 
at oo; that is, for c > 0, j{x)~^dx = oo. When there is killing at 0, the process is recurrent 
only when the scale function is infinite at both ends. In analytic terms, recurrence means 
that the associated generator is critical (see [19] and [10]). Recall that L'^ is called critical iff 
there exists a unique (up to constant multiples) positive solution -0 of L'^ip = 0. Otherwise 
L'^ is called subcritical. We know from criticality theory — for example, from Theorem 3.15 
of [19] — that the generator must be critical if is an isolated eigenvalue. A generalization 
of this fact will be used in Lemma 13.31 

The semigroup e~*^" has a probabilistic representation: We consider the product space 

C7([0,oo)) X [0,oo) = {{uj,0 G C([0,oo)) x [0,oo)} 

endowed with the natural product cj-field. Let (PxOxe(o,oo) denote the family of measures 
which is induced by the Dirichlet form q^. For x G (0, oo) we define the measures 

Px ® e"^ d^. 

and the stopping time 



r,(6^,0 =mf|s > I ^ K(a;,)(is > e|. 



If we set 

TQ = min(ro,r«;) 
then we have the Feynman-Kac representation 

(15) (e"*^V)(a;) = [fiXt),T9 >t]=E, [e-/o «(^^)'^7(^t),To > t 
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Since we are working with the self-adjoint version of the generator (with respect to the measure 
r) , the Feynman-Kac representation holds in great generality, following the derivation in [16] . 
We will generally omit the tilde, since it will be clear from context which measure is meant. 

Let us recall the usual Feller classification (see e.g. chapter 3 in [32]) of boundary points 
for diffusion generators — b{x)-^ in an open interval (0, r). 

Definition 2.3. Let c G (0, r) be given and set 7(x) = e/"2%)'ij/. The point r is called 
accessible, if ^(y) j{z)^^ dz dy < oo, and otherwise inaccessible. If r is an accessible 
boundary point, then it is called regular iff j[x)~^^{y) dx dy < oo. If r is accessible 
and jiu) {f^ dx^ dy = oo, then r is called an exit boundary. If r is inaccessible 

then it is an entrance boundary, iff ^{x)^^^{y) dx dy < oo. If r is inaccessible and 
jl l{v) iJc "fi^)'^ ~ then r is called natural. Of course the same classification 

holds for 0. 

Except where otherwise indicated, we will always assume that the boundary point oo is 
inaccessible. 

It is easy to check that the boundary point r is regular if and only if 7(2;) dx < 00 and 
'y{x)~^ dx < 00. A boundary point is thus regular in the sense of Feller if and only if it 
is regular in the sense of Weyl (c/. [21] )• Let us recall the relevant definition from the Weyl 
theory of selfadjoint extensions of singular Sturm-Liouville operators L'^ = —^^(7^) + 
in (0, r), adapted to our special situation. 

Definition 2.4. Let z £ C We say that boundary r is of limit-point type, if there exists 
c E (0,r) and 2; G C, and a solution f of {L'^ — z)f = such that \f{y)\'^'yiy)dy = 00. // 
there exists c G (0, 00), such that for every solution of the equation (L'^ — z)f = the integral 
|/(y)p7(y) (iy is finite, then we say that r is of limit- circle type. An analogous notation 
applies to the boundary point 0. 

A fundamental result in the theory of Sturm-Liouville operators is the so called Weyl- 
alternative, which states that exactly one of the above situations holds and that the limit- 
point/limit-circle classification is independent of z G C (see [M])- Moreover if we are in the 
limit-point case at r then for every C \ M there exists exactly one solution of the equation 
{L'^ — z)f = which satisfies /J |/(s)p 'y{y)dy < 00. Roughly limit-circle case at a boundary 
point r means that we have to specify boundary conditions at r in order to get a selfadjoint 
realization, whereas in the limit-point case at r no boundary conditions at r are necessary. 



2.3. Quasi-limiting and quasi-stationary behavior. We say that Xt converges from 
tlie initial distribution to the quasistationary distribution ip on compacta if for 

any positive z, and any Borel A C [0,z] 

t^oo (p{y)-f{y)dy 

Xt converges from the initial distribution v to the quasistationary distribution ip 

if /q°° ^p{y) "f{y)dy, and for any Borel subset A C [0, 00) 

hm F^{Xt £A\Ta>t)= 00 . ^ . ^ , ■ 
t^oo J (f{y)j{y)dy 
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Finally we say that escapes from the initial distribution i/ to infinity if 

lim F^(Xt <z\Tg>t)=0 

Remark 2.5. In the literature there is no completely standard terminology for quasistationary 

distributions. The probability ^ ,{ , described here is sometimes also called a 

Jo fiyjiiym 

quasi-limiting distribution. A quasistationary distribution v is often defined as a probability 
measure v supported in (0, oo) satisfying 

Fi){Xt £ A\T9>t) = i>{A), V Borel sets A C (0, oo), t > 0. 

Quasilimiting distributions are also called Yaglom limits. It is not difficult to see that quasilim- 
iting distributions are also quasistationary distributions. 

2.4. Previous results. Observe that we have in equation ([7]) that (/'(-^O) ■) is positive. Stein- 
saltz and Evans ^47j showed 

Theorem 2.6 (Theorem 3.3 in [17]). Assume that oo is a natural boundary point and that 
we are in the limit-point case at oo. Suppose that either 

liminf k;(x) > Aq or limsupK(a;) < Ag. 

X—KX) X— !>00 

Then either Xf converges to the quasistationary distribution '/'(Aq, •), or Xf escapes to infinity. 
In the case liminf k(x) > Ag, Xf converges to the quasistationary distribution (/^(Ag,-) if and 
only if /g°° (p{Xq, y) -i{y)dy is finite. 

A priori it would not have been clear that the conditional distribution converges, and 
that the mass cannot split, with part of the mass remaining on a compact interval and the 
remainder escaping to infinity. Having recognized that there is a a dichotomy, it is natural 
to then seek a simple criterion for discriminating between the cases: escape or convergence. 
One such is given in [37], under which Xt converges to quasistationarity, namely when Ag < 
K =: limt_s.oo '^(i) together with the growth bound 

{GB') 36, K > Vy large enough : \b{y)\ < by and K{y) < ky 

or the related bound 

{GB") 3bi, 62, K, /? > Vy large enough : 61/ > b{y) > -biy, b'{y) > -hy"^ 

and K(y) < ky 

While these conditions are satisfied in many applications they are from a theoretical point of 
view unsatisfactory. In particular, it seems peculiar that an upper bound on the killing rate 
as in {GB') should be necessary. On the contrary increasing the killing rate /t should, from a 
heuristic point of view, only strengthen the convergence to quasistationarity. 

Remark 2.7. We make use of Theorem \2. 6\ only in the case Ag > lima;^oo ^(x) andT{{0, 00)) < 
00. In the other cases we use different techniques. In the next chapter we will show that 00 is 
always in the limit-point case. As emphasized in [47J in this case the heuristic behind Theorem 
\2.6\ is quite clear, but the translation of this idea into formal mathematics is not trivial. 
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3. Analytic Results 

In this chapter we derive several key analytic facts about the spectra of generators and 
resolvents. While some of these are standard in the theory of Sturm-Liouville operators, and 
well known to specialists in that field, they are less familiar to probabilists, and we explain 
them in some detail here. In section [3?T] we show that the technical conditions for a limit-point 
boundary at oo may be weakened. Section 13.21 derives basic results linking the spectrum and 
speed measure. Section 13.31 presents the standard parabolic Harnack inequality in the form 
that we will be using. Section 13.41 applies the analytic results to convergence on compacta. 
Section [3.51 explains why strong conditions on the initial conditions are unnecessary. Finally, 
section \3M generalises the results to the case of an entrance boundary at cxd. 

3.1. Classification of boundary points. We start by establishing a connection between 
the Feller classification and the Weyl classification of boundary points. This has already been 
investigated in |55j for the case k = 0, but in this work the author introduces the notion of 
weak entrance boundary and shows that one is in the limit-circle case if the boundary point is 
of weak entrance type. We show that there are no weak entrance boundaries at oo by proving 
that oo is in the limit-point case if and only if it is not a regular boundary point. The proof 
we give is well known in the Schrodinger case (see ^ for similar ideas in a much more general 
context). We assume regularity of the coefficients of the Sturm-Liouville expression, although 
weaker assumptions would also suffice. 

Lemma 3.1. Let the Sturm-Liouville expression Tf{x) = —^^i^{j{x)f'{x)y + K{x)f{x) he 

given. Assume that 7 is strictly positive and locally Lipschitz in (0, 00) and k G £^^^^([0, 00)) 
such that k{x) > — C|xp + D for some constants C,D > 0. Then we are in the limit-point 
case at 00. 

Proof. We can assume, without loss of generality, that D = and that 7 is continuous up to 
the boundary. The first assumption is obviously harmless. If 7 is not continuous up to zero 
we can consider the differential expression in (l,oo) instead of (0,oo). This shift does not 
change the Weyl-classification of r at infinity. Similarly, we may assume that the boundary 
condition at is Dirichlet {a = 00), since the classification at infinity is unaffected by the 
boundary condition at 0. 

As usual in the theory of Sturm-Liouville operators we define the maximal operator T and 
the minimal operator T associated to the differential expression r as 

P(T) := {/ G £^ I /, 7/' absolutely continuous in (0, 00), r/ G 

Tf := Tf for / G V{T) 

and 

V{f) := { / G V{T) I / has compact support in (0,oo)}, f / : = r/ for / G V{f), 
respectively. Let To be the restriction of the maximal operator T to the domain 

p = {/GP(r)|/(o) = o}, 

i.e., we put Dirichlet boundary conditions at the boundary point 0. 

The deficiency indices of Tq are (1,1) if the limit-point case holds at 00, and (2,2) if 
limit-circle holds at 00. In the former case, the maximal symmetric extensions of Tq are one- 
dimensional; in the latter case, they are two-dimensional. If defines a symmetric operator 
— (/, To/) G M for every f £ V — then it cannot be limit-circle, since To would be a 
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symmetric three-dimensional extension of Tg, whereas the maximal symmetric extensions are 
two-dimensional. 

It remains to show that Td is symmetric. Let (p £ C^(M) such that < (/? < 1 and 



ip{x) 



1 if |x| < 1 
if|x|>2. 



Further we set (fkix) = (^(|) {k G N). This gives, for / G P(rD) and k gN 

POD 

(/,rz,/)=lim / ipk{xff{x)TDf{x)-f{x)dx 
k^'^ Jo 

' 1 

Mxffix) [-^ilf'nx) + k(x)/(x)] 7(x) dx 



(16) 



hm 
hm 
hm 

fe— )-oo 



ipk{xffix)y-fix)dx+ / ipkixfn{x)\f{x)\'^'y{x)dx 



dix)Uf'ix)\' + K{x)\f{x)\')^{x)dx 



+ 



ipk{x)(p'f,{x)f{x)f'{x)'y{x) dx }. 



Observe that in the third hne the boundary term iflfjf |o° coming from the integration by 
parts vanishes, since 7/ is continuous up to (since it satisfies an ODE), /(O) = and ^Pkix) 
is identicahy zero for x large enough. 

The first term on the right-hand side is real and we have to prove that the second term 
converges to as A; — )• 00. We have, by the Cauchy-Schwarz inequality and the properties of 
the cut-off sequence {(pk) 



(pk{x)ip'k{x)f{x)f'{x)-f{x)dx 



(17) 



< 



< Ck 



•fkixflf 



/■c 

' {x)\'^j{x) dx 

Jo 



\ip'kix)f\f{x)\'^{x)dx 



(/•oo /•2k 
/ Mx?\f'ix)\^^{x)dx \fix)\^^{x)dx 
Jo Jk 



For the first integral on the right-hand side we integrate by parts in a similar vein to (|16p . 
The assumptions on k as well as the elementary inequality \ah\ < /A + lP' imply 



^k{xf\nx)\'^{x)dx 



^k{x) f{x){TDf{x))j{x)dx 



ipkixfK{x)\f{x)\'^-f{x)dx 



^k^'k{x)f{x)f'{x) 'j{x)dx 



< 



Mx)\f{x)\\Tnfix)\-f{x)dx + C I iplx^\fix)\-'-fix)dx 



+ 



\ipk{x)ip'k{x)f{x)f'{x)\ -f{x)dx 
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< 



< 



Mx)\f{x)\\TDf{x)\j{x)dx + C / ^ixV{x)\'^{x)dx 



+ 



-\Mx)nxr + Mx)fixr 



'y{x)dx 



\\TDf\\ + C{2kf 

^k{xf\nx)\'j{x)dx + M'k-^ 



1 

+ 4 



This yields 



Mx?\nx)\' j{x)dx < i\\f\\\\Tgf\\ + 16CA;2||/||2 + 4M^k-' 
and therefore for large k 



(18) 



ipk{xf\f'{x)\'' dx<Ci + C2k' < Csk'. 







Thus the inequalities pT|) and ([T8]) impliy that (observe that / G L^((0, oo), 7)) 

r2k 



ipk{x)ip'k{x)f{x)f'{x) j{x)dx 



< Ck-^ f Csk"^ 



|/(x)| -f{x)dx 



0. 



as A; —7- 00. This proves the assertion, and so completes the proof. 



□ 



3.2. The spectrum of Sturm-Liouville operators. We begin with a version of the spec- 
tral theorem for self-adjoint operators on a Hilbert space, specifically adapted to Sturm- 
Liouville operators. A proof of it can be found in general references on the theory of Sturm- 
Liouville or Schrodinger operators, such as pOj [HI |58]. 



Let r 



j_ d 

'27 



^(7^) + K be a Sturm-Liouville expression which is regular at and in the 
limit-point case at infinity, and let H be the self-adjoint realization of r in £^ with boundary 
conditions [TH at 0. Let (p{z, •) be the unique solution of the ordinary differential equation 
Tip{z, •) = zip{z, •) satisfying (p{z, 0) = 1/(1 + a) and ip'{z, 0) = a/(l + a). 

Given a continuous function F G C(M), we have a corresponding maximal multiplication 
operator Mp on i2^(M, ^) defined by 

V{Mf) = {g£ £'(M,//) s.t. gF G ^'(IR,^)}, 
Mpig) = Fg. 

Theorem 3.2 (Weyl's Spectral Theorem). There exists a measure a whose support is Ti(H), 
such that the map taking a compactly supported function h 2? to the function h G 2? {Tj{H) , a) , 
defined by 

/•oo 

^(■) = / h{x)(p{-,x) 'j{x)dx 
Jo 

may be uniquely extended to a unitary mapping U : £2((o,oo),r) £2 ^/jg property 

UF{H)U'^ = Mp. 
The spectrum of H is simple, and S(F(iJ)) = essrana{F). 
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The spectrum of a self-adjoint operator A may be divided into two components: the 

essential spectrum 'Eessi-A), comprising the limit points and eigenvalues of infinite multiplicity; 
and the discrete part E^(A), comprising the isolated eigenvalues of finite multiplicity. In the 
Sturm-Liouville case every eigenvalue has finite multiplicity (no more than 2), so the essential 
spectrum consists only of limit points of the spectrum. It is well known that the essential 
part of the spectrum of self-adjoint operators is invariant with respect to relatively compact 
perturbations (see Theorem 9.15 of t53j). (We recall that an operator V : X ^ X on the 
Banach space X is called relatively compact with respect to T : X — )■ X if T){T) C ViV), 
and if for some z G C \ S(T) the operator V{T — z)'"^ is compact. We refer to section 9.2 of 
|53] for further details.) 

The core of our results is contained in the following analytic lemma, which catalogs some of 
the key linkages among the base of the spectrum, the scale measure, and the speed measure. 
These take us beyond the results of Theorem l2.61 by separating the influence of the drift from 
the effect of the killing term. Moreover, they show clearly why the case \q < K will turn out 
to be easier than the case Aq > K. The major results — particularly Theorems 13.151 [4.11 14-71 
and l4.9l — will in essence be just unpacking these analytic results in probabilistic terminology. 

Lemma 3.3. With the above definitions, 

(i) //lim^^oo k(x) = K, then Se.s(i°'") +K = Se.s^L'^'"). 

(ii) Ao > and j{x)~^ dx = oo imply r(IR+) = ^{x)dx < oo. 

(iii) Ao > and r([0, oo)) = oo imply Xq'^ > 0. 

(iv) Ao > and r([0, oo)) = oo imply lim,.^oo f logr([0, r)) > 0. 

(v) // Aq < liminfj^^oo k{x) then Aq is a simple isolated eigenvalue with a unique positive 
eigenfunction. 

(vi) If a > (not pure reflection at 0) or if T is infinite, then is not an isolated 
eigenvalue of L^'". 

(vii) If a > (not pure reflection at 0) or if T is infinite, then Aq > limsup^_^cxD '^(^) 
implies Aq" > 0. 

Proof. Assertion can be derived from the fact that the essential spectra of two self-adjoint 
operators Li and L2 coincide if for some z £ C \ (E(Li) U £(^2)) the difference 



is a compact operator (c/. Theorem 9.15 of [53j)- Set Knit) = 1[q^„] (t)(K(t) — K). The 
resolvent equation gives for 2; G C \ M 



(L"" - z)-'^ - (L - zy^ = (L"" - z)~^(L - L^"){L - 2)"^ = -{L^" - z^^K^iL - z)-^. 



Observe now that the operator k„(L — z)~^ is compact; that is, the operator acting by 
multiplication with k„ is relatively compact with respect to the operator L. This can be seen 
by considering the explicit form of the resolvent (see chapter 3.3 in [21]; similar results can 
be found in [II]). We have 



(Li - z) 



1 



{L2 - z) 



1 



[k„(L - z) '^]g{x) 



Kn{x) 



1 



v{x) I u{y)g{y)j{y)dy + u{x) / v(y)giy)j{y)dy 



W{v,u) 
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where u and v are linearly independent solutions of 
(r — z)w = satisfying 

«(0) = u{0) = — and 

1 + Q 1 + Q 

b(y)p7(y)f^y < oo- 



Observe that here we use the fact that we are in the limit-point case at infinity. The Wronskian 
W{f,g) of two locally absolutely continuous functions / and g is defined by 

W{f,g){x) = [f'{x)g'{x) - f'{x)g{x)] j{x), 

Thus Kn{L — z)~^ is an integral operator in 2,'^ with kernel k{-, •) given by 

W{v,u)~^Kn{x)v{x)u{y) if y < X, and 
W{v,u)~^Kn{x)v{y)u{x) if y > X. 

The known properties of u and v imply 

\k{x,y)\'^ 'y{y)l{x)dy dx 



k{x,y) 



JO 



\uiy)\'^ l{y)dy + j \v{y)\^ l{y)dy\ \Kn{x)\^-i{x)d. 



Thus — z)~^ is Hilbert-Schmidt, hence also compact. 

We complete the proof by observing that the resolvent equation 

implies 

IKL-'^ - Z)-^ - (L- - < - z)-^\\\\k -Kn- i^lloolKL" - z)-^\\ 

[xszY 

as n — )• oo; that is, L'^"~^^ converges in the norm-resolvent sense to L'^. In the second in- 
equality we used the fact that the operator norm of the operator that acts as multiplication 
by a function / is just the supremum norm of /. 



X < oo. 



Assertion (In]) is contained in [35] and also follows from Theorem 1 of the recent work [IT] . 
(In ^41j somewhat stronger conditions on the drift are imposed, but these are not actually 
necessary for the proof.) 



Assertion (imh : Because L^'°^ and L^'^ differ only in their (one-dimensional) boundary con- 
ditions, the difference (L°'° + 1)-^ - (L°'° + 1)"^ has one-dimensional range, so is compact. 
(For more details, see [53^ Satz 10.17].) The bottom of the essential spectrum of L^'^ is then 
strictly positive, since it coincides with the bottom of the essential spectrum of L^''^, hence 
is above the bottom of the full spectrum of L*^'". If Ag'" := inf spec(L'^''') = then Aq'' = 
is necessarily an isolated eigenvalue of the operator L^'^. Let us assume that Aq*^ = 0. The 
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unique (up to positive multiples) non-trivial and non-negative eigenfunction V]\f G £ associ 
ated to Aq*^ = therefore solves the boundary value problem 



L''%j, = A°'%^ = 0, vj,{0) = ^ and 7^(0) = 



a 



l + a 



Since this ordinary differential equation has a unique solution, and since the constant function 
1 is also a solution of this equation, we conclude that Viv = 1. Thus 1 £ 2p, which means 
that r((0,oo)) < 00, contradicting our assumption that T is infinite. It follows that Ag' > 0. 



Assertion (|iv|) follows from the above and the work of Notarantonio [3H]. His result implies 
that the bottom of the essential spectrum of the operator L^'^ with Neumann boundary con- 
ditions at is bounded above by lim sup^_j.(,o ^ logr((0, r)). This is if the volume growth is 
subexponential. Since we have already showed that Aq'^ > 0, the result follows. 



Assertion (jvj): Assume first that lim^^^oo ^^(2;) exists. If Aq < lim2._j.00 = K then an 
application of the result @ shows that L'^ = L + K + (k — K) has the same essential spectrum 
as L + K. Since L is a positive operator, the bottom of the spectrum of L -\- K, hence a 
fortiori of the essential spectrum, has to be at least K, hence bigger than Aq, which implies 

m- 

Let us now assume only that liminf^^^oo /^C^;) > Aq. By the decomposition principle (see 
section 131 in [3]) it is not difficult to see that L'^ has the same essential spectrum as the 
operator (a > 0), defined as the selfadjoint extension of t'^ in £^((0, oo),r) satisfying 
Dirichlet boundary conditions at a. If oq > and e > are such that 'mix>ao 't(a^) > Aq + e 
we conclude that 

infSess(i^") >mfS(L^J 

> inf / \(f' {x)\'^ ^{x)dx + / K{x)\ip{x)\'^j{x)dx 

ll'^ll£2((a(,,oo),r)=l 

>AS + e. 



Assertion (jvij): Let f,g be continuous functions, nonnegative and not identically zero, with 
compact support on (0, 00), and Ai := inf S(L'''")\{0}, where a > 0. Suppose is an isolated 
eigenvalue, so that Ai > 0. Then 

(f,e-'^%)= [ e-'^d{f,Eg)iX) 

/•OD 

= {f,Em)g)+ / e~'^d{f,Eg){X) 

{f,uo){9,uo) as t 00. 
This is positive, since uq may be chosen to be strictly positive. 
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By Lemma |2 . 1 1 there is a constant C such that for all x G supp(/), 
|e-*^°5(x)| <C||VZ0e-*^"5|| 

1/2 



so that 



< C||(7||e~'^^* for t sufficiently large, 
/, e-*^^) < C\\g\\e-^'' f{xh{x)dx ^ 



as t — 7- oo, which is a contradiction. 



Assertion (jviip : Suppose first that the limit K of k(x) exists. Intuitively, what we are 
saying is that when the mass in a neighborhood of shrinks at a rate faster than K (what Aq 
measures), it is being driven by drift: Either the mass is being swept down into a region of 
high killing near 0, or it is being swept up away from 0. In the latter case, the drift will still 
cause the mass near to shrink exponentially in the absence of killing; in the former case, 
the killing at will do the job, except in the case of pure reflection at 0. 

By part (ji]), we see that = L + K + (k, — K) and L + K have the same essential 
spectrum. In particular we conclude that inf Se<js(-L) + K = inf Sess(-L + K) > Ag and 
therefore inf $]ess(L) > Aq — > 0, so Aq = would imply that is an isolated eigenvalue. 
Since this is impossible, by assertion fvH) . it follows that Aq > 0. The extension to the case 
when the limit does not exist goes exactly the same way as in the proof of assertion (jvj) 
above. □ 

Remark 3.4. R. Pinsky has showed |41j that the conclusion ^ of Lemma 1 3. 3\ can be sharp- 
ened. Assuming that absorption is certain, it was shown that 

(19) ^TTTIT < ^0 < 



8A{b) - - 2A{b) ' 
where 

A{b) = ( r ^{V)dy\ ( [%{y)-'dy 



Related analytic inequalities, which are usually referred to as weighted Hardy inequalities, can 
be found in [37\. Indeed the results of [37J imply Pinsky 's bounds. 

Remark 3.5. The fact that the bottom of the spectrum is an isolated eigenvalue is also of 
practical interest, because in this case the associated eigenfunction can be approximated ac- 
curately by the ground states of regular Sturm- Liouville operators on bounded intervals (see 
the recent survey [54] j. Such a result has recently been rederived [51\, in the context of ap- 
proximating the minimal quasistationary distribution of a diffusion generator with discrete 
spectrum via interacting particle systems of Fleming- Viot type. 

Remark 3.6. Assume that Ag is an eigenvalue with associated eigenfunction ux^ G which 
by general theory is strictly positive and simple. Then 

(20) lim e^otp^{t,x,y) = cux-{x)ux^{y), 
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where c is a normalizing constant. This was proved in |45| for the transition function of Brow- 
nian motion on Riemannian manifolds hut the proof carries over without essential changes to 
our case. 

We will also make use of the following result which is a special case of Theorem 3.1 in [47]. 

Lemma 3.7 (Theorem 3.1 in [47j). Let < f £ 2,^ with compact support supp{f) C [0, oo) be 
given and let uj denote the measure f{x)^{x)dx. Let L'^ be as in LemmalKM and let p'^(t, •, •) 
denote the integral of e~^^'^ . Then for arbitrary measurable bounded sets A,Bc (0,oo) 

^.^ /o°° fjx) Ib P^i*^ ^' y) ^(y) dy dx ^ f^(p{X'^,y)j{y)dy 

f(^) lAP'^(*^^^y)^(yh(^) dy !A'^iK^y)iiy)'^y' 

i. e. Xt converges from the initial distributions jao ^^^^^ on compacta to the quasistationary 
distribution <y9(AQ, •). 

The above Lemma can be proved directly, using the spectral representation for Sturm- 
Liouville operators. The reader will see the necessary arguments later in this work in the 
proof of Theorem 13.151 Our first goal is to extend this result to the case of general compactly 
supported initial distributions v. 

We begin by deducing some consequences of Lemma 13.71 This will lead to Proposition 
13.91 which is a "strong ratio limit theorem". Before we start proving the strong ratio limit 
theorem we explain another analytic fact which has no direct relation to spectral theory but 
which will turn out to be very useful. 



3.3. Local parabolic Harnack inequality. A crucial tool for smoothing analytic infor- 
mation about the transition kernel between different times and sites is the local parabolic 
Harnack inequality, which quite general holds for second order parabolic differential equa- 
tions. One version appropriate to our current purposes may be found in [3l], and states that 
for fixed xo,i G (0,oo) and i? > there is a constant C such that for every weak solution u 
of [dt + L'^)u = which is non-negative in Q{{xq^ to),AR) C (0, oo) x (0, oo) 

sup u < C inf u, 

e{(xo,to),R/2) Q({xo,to),R) 

where 

Q{{xo,to), R) = {X £ R |max(|a; — xo \,V\t-to\) < R,t< to} 

and @{{xQ,to), R/2) = Q{{xo,to — R'^),R). As in Theorem 10 of [Hj this inequality can be 
applied to the transition kernel p'^{t, x, y) in order to prove that for every compact K C (0, oo) 
and T > there is a constant c = c{K, T) > such that for t >T, xi,X2,X3,X4 £ K 

(21) C-^p'^it, Xi,X2) <p''{t,X3,Xi) < Cp'^{t,Xi,X2). 

Moreover the local parabolic Harnack inequality shows that there exists a locally bounded 
function ( : (0, oo) — )• (0, oo) such that for every t > 1, y > 0, and x,z > satisfying 

|z — xl < ^ A ^ 



(22) 



p''{t,x,y) < C{x)p''it + l,z,y). 
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3.4. Strong ratio limit theorem and convergence on compacta. 
Lemma 3.8. For any fixed xq £ (0,oo) the family of functions 

[0,oo) X M+ X M+ 9 {t,x,y) ^ ^"f^^"'"^'^^ | s > 1 

is relatively compact in the space C((0,oo)^,M) of real-valued continuous functions on (0, cxd)^, 
endowed with the vague topology. 

Proof. Let {sn)neN be a sequence with 1 < s„ — )• oo and set for t G [0, oo), x,y £ (0,oo) 

X P'^{t + Sn,x,y) 
rn{t,x,y) = — — -, 

p'^{Sn,Xo,Xo) 

where a £ (0,oo) is fixed. The functions (t,x,y) i— )■ r„(t, (n G N) are solutions to the 
parabolic equation 

(29t + LS + L^)r„(t,x,y) = 0, 
where the operator and Ly act as L'^ on x- and y-variable, respectively. By the local 
parabolic Harnack inequality (see section 13. 3p we conclude that for each compact set K C 
(0, oo) there exists a constant Ck such that for all n G N, t > and x,y,a G K 

p'^{t + Sn, X, y) < Cxp'^it + Sn, Xo,Xo) 

By general spectral theory it is proved in [H] that r i— )■ p'^{r,XQ,XQ) is non-increasing. There- 
fore we conclude that for t > and x,y £ K 

p''{t + Sn,x,y) 

p''(Sn,Xo,Xo) 

Theorem 6.28 in [31j shows that the set {r„ | n G N} is locally uniformly equicontinu- 
ous. Therefore by the theorem of Arzela-Ascoli |26t Theorem 17] there exists a subsequence 
{rnf,)k&i which converges locally uniformly. □ 

This proof is modeled on Theorem 2.2 of [3]. Since that theorem assumed the operator was 
critical and the coefficients were Holder-continuous, some modification was required 

The analytic core of quasilimiting behavior is the convergence of ratios of transition kernels, 
which we state and prove here as Proposition 13.91 This will imply convergence to the qua- 
sistationary distribution on compacta. Theorem 13.151 Convergence on the whole state space 
will then require a consideration of the recurrence or transience, to decide whether most of 
the mass stays in a compact interval or escapes to infinity. 

Results comparing transition probabilities at different times and sites, in the limit as time 
goes to infinity, are commonly referred to as strong ratio limit theorems. Strong ratio limit 
theorems for certain branching processes can be found in [5]. A proof of the strong ratio 
property for certain Markov chains on the integers was given in [27j . 

Proposition 3.9. For any a G (0, oo) 

^.^ p{t + s,x,y) ^ ^-A«t V^(Ag,a;)v9(Ag,j/) 
^*^oo p{s,a,a) (/j(Ag,a)v3(Af},a)' 

Proof. For every sequence (sn)neN C (0, oo) converging to infinity we know by Lemma 
that for some subsequence {snk)k of (sn) there exists a function ip such that 

p'^it + Sn^.x.y) 



p'^{sn,^,a,a) 



ip{t,x,y), 
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where the convergence is locally uniform in [0, oo) x (0, oo)^. Since by Lemma 7.5 in [37] (see 
also [H]) for every / S -G^ with compact support 



one easily concludes that 

ij{t,x,y) = e-^o(.i^{0,x,y). 
Lemma 13.71 shows that for every f,g,h £ C^(0, oo) 

fi^ 9{y)v^iK^y)^iy)dy ^ ^.^^ fi^) /o°° 9{y)p''{snk , x, y) i{y)i{x) dy dx 

/o°° h{y)'P{K^y) l{y)dy fix) h{y)p'^{snk > x, y) l{y)l{x) dy dx 

_ ^. /o°° /(^) /o°° g(?/) p^(f!^txo,'x!) 7(y)7(^) dy dx 

" '^^^ r /(-) r 7(y)7(-) 

_ Jo°° /(^) /o°° 9{y)ip{0, X, y) l{y)l{x) dy dx 
/o°° /(^) /o°° ^(y)^(O' 2;, y) 7(^)7(2:^) dy dx ■ 
This implies that for x G (0,oo), g,hG C(f ((0, 00) 

9{y)v{>'o,y)i{y)dy f°° f°° 

c^ui \ t\K \ / \^ / Ky)'>P(.0,x,y)-/{y)dy = g{y)ilj{0,x,y)j{y)dy, 

Jo KyMK^y)i{y)dy Jo Jo 

and hence for every h 

Vy^.x.y) = ip{Xo,y) 00,. ^ ^ . X , 

Due to the symmetry oi ip{0,-,-) we conclude that for some constant c > 

ij{0,x,y) = cv?(A;j,x)v?(AfJ,y). 

Because of ^{0,a,a) = 1 we arrive at = (^(Aq , a)(^(AQ , a). Since this is true for every 
subsequence, the assertion of the theorem is proved. □ 

Corollary 3.10. // 1/ is any compactly supported initial distribution, and f a nonnegative 
compactly supported measurable function with > 0, then for any fixed t, 

|logE,[/(Xt+,)]-logE,[/(X,)]| 

is bounded for s G . 

Proof. Since g{s) '■= J J p'^{s, x, y)'^{y)f{y)dv{x) is positive and continuous, it suffices to show 
that 

n/i- • , 9{t + s) g{t + s) 
< limmt — — < limsup -j-^ — < 00. 

g{s) s^oo g{sj 



By ([2T]) we may find positive c such that for s sufficiently large, and any x*,y* G K 
supp(i/) U supp(/), 

c-V(s,x*,y*)r[/] < g{s) < cp^{s,x^,y,)r[f]. 

Thus, 

^_2 P'^{s + t,x^,y^) ^ g{s + t) ^ ^2 P^is + t,x^,y^) 
p'^{s,x^,y^) ~ g{s) ~ p'^{s,x»,y^) 
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By Proposition 13.91 the upper and lower bounds converge to c^e~'^o and c~^e~'^o respectively, 
as s — )• oo. □ 

Remark 3.11. In terms of parabolic Martin boundary theory Proposition \3. 9\ says that every 
sequence {sn,x) C (0, oo) x (0, oo) with lim„_s.oo Sn = oo converges in the parabolic Martin 
topology to the parabolic Martin boundary point corresponding to the minimal parabolic func- 
tion hx^{t,x) = e^o^ ip{Xq , x) . The parabolic function hx^ must actually be invariant, since it 
corresponds to a point in the parabolic Martin boundary whose time coordinate is oo. 

Remark 3.12. The existence of strong ratio limits for general symmetric diffusion, i.e. the 
existence of 

r P{t,x,y) 
hm — -, 

t^co p{t,xo,yo) 

where p{t, •, •) denotes the transition kernel of the diffusion, was investigated under special 
conditions by Brian Davies in [14J and is now often referred to as Davies' conjecture. In a 
private communication Gady Kozma disproved this conjecture by presenting a counterexample. 
Proposition Ig.gl shows that in one dimension the Davies conjecture is true, if one boundary 
point is regular. It is an open question, whether the Davies conjecture generally holds in one 
dimension. 

Proposition 3.13. For all positive z, including z = oo, 

lim sup - log Xt < z\ and lim inf - log F^^Xt < z\ 

t^oo t ^ t^oo t '- ^ 

are both constant in x > 0. Hence also 

limsup — logP^jXt < z\tq > t\ and liminf — logP^j^t ^ zItq > t\ 

t^oo t ' ^ t^oo t ' ^ 

are both constant in x > 0. 

Proof. We prove only the first statement for limsup, the other proof being identical. Suppose 
we have x, x' such that 

(23) limsup - logP^.{Xt < z} < limsup - logF^>{Xt < z). 

t— ^oo t t—^oo t 

We may assume without loss of generality that \x — x'\ < ^ A A (If not, then there 
must be other starting points closer together where the limits differ.) Applying ()22p . we have 
for all t > 1, 

P.'{Xt<z} <Cix')F4Xt+i<z}, 

so that 

1 r T 1 r 1 1 F^{Xt+i < z] 

limsup — logPj^/jXj < zj < limsup — logPx-j-'^t < ■zj + limsup — log ■ 



t—^oo 



t t^oo t t~^oo t rx{Xt < Zj 



The limits on the second line are by Corrolarv 13.101 We have then a contradiction to (j23p . 
which completes the proof. □ 
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3.5. Conditions on the initial distribution. In their version of the ratio hmit theorem 
\47\ Theorem 3.1], Steinsaltz and Evans had to pose an additional condition on the initial 
distribution and they stated the general case as an open problem. Their most general 
condition reads 

{ID') If Xq has distribution v, then 3s > for which the distribution of Xg 
has a density / G with liminf f7/(A) > -co. 

It is not obvious how to verify that a given initial distribution with compact support satisfies 
the condition {ID'). Using some results from spectral theory and several ideas of [47j, we can 
remove this restriction. An essential ingredient in the proof is Lemma l2.H which allows us to 
ignore the upper end of the spectrum for large t. 

Lemma 3.14. Given g G 2,'^ , we have 

(24) lim t"^ log ||e-*^''5|| = -A.. 

t— )-oo 

Proof. By the spectral theorem ([5]) we know that 

2 

limsupfMog lle"*'^"^!! = limsupfMog / e^'^^^dWE^gW {X) 



^9 

< — 2Ag + limsupt log ll^lf 
= -2A„ 



'9 

For the lower bound we take any A* > Ag, and have 

r-A. 



liminft-Mog||e-*^"c,||^ > liminft-^log / * e'^^^dWE^gf {X) 

t—>-00 " " t—KX) 

> -2A, +liminft-ilog||^''([Ag,A,])5f 



Since Xg is in the support of this is equal to — 2A*. Since this is true for any A* > A^, 

this completes the proof of (p^ . □ 



For a Radon measure on (0, oo) and a Borel measurable function / : (0, oo) i— )• C we use 
the notation (i^, /) := f{s) ^{ds). 

Theorem 3.15. The killed diffusion Xt converges on conipacta to the quasistationary distri- 
bution with density proportional to ^p{Xq, ■) from any initial distribution which is compactly 
supported in (0,oo). 

Proof. An application of Weyl's eigenfunction expansion theorem and Fubini's theorem tells 
us that the operator ^^'^([Aq, Ai))e~*^"' has a continuous integral kernel 

(25) {t,x,y)^h^'{t,x,y)= [ e-'^y^{X,x)^{X,y) da{X) 

with respect to the measure F. This implies that for every compact subset K C [0, oo) the 
function E'^{[Xq, Xi])e~^^'^ Ik is continuous and therefore 

(z.,S'^([AS,Ai])e-*^^li^)= [ E-{[X^,Xi])e-''^\K{x)du{x) 
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is well defined. For every Borel set A C [0, z], then, 



(26) 



-t\ 







Let g : [Ag,oo) R 
(27) lim sup 



[AS,Ai] 

be any continuous function. Then 



ip{X,x)di'{x) / ip{X,y)dT{y) 



da{X). 



t—^oo 



/[As,A.]e-*"5(A)da(A) 



da{X) 



< sup \giX) - giX^)\. 

[A^.Ai] 



As in the proof of Theorem 3.1 of [U], for any Ai, Ai, A2 > Ag set A* = Ai A Ai A A2 A A2 and 
A* = Ai V Ai V A2 V A2. Then we have the bound 



/[AS,A,] ^-''9{X) da{X) .^^j e-^'g{X) da{X) 



'[AS,A2] 

which tells us that 



da{X) 



(28) 



lim sup 



[A^.Ai] 



^AS,A2] ' 

r''^g{X)da{X) 



-t\ 



da{X) 



,H/[ Ag,A>] lg(A)|rf^(A) 

/fAS,A.)l MX) 



Since g is continuous, ()27p and (j28p combine to show that 



-tx 



da{X) 



is independent of Ai,A2 € (Ao,oo). 



(29) 



lim 



/[Ag,A,]e"*'5(A)da(A) 



[AS,Ai 



da{X) 



By (i26]l we now see that for every Ai € (Ag, 00) 
iy[E%[X^,Xi])e-''^^lA] 



lim 

t—>-oo 



(30) 



i.[i?«([Ag,Ai])e-*in[o,,]] 

/[As,Ai] e"*^ [/o' ^) Ia <^(-^' y) ^r(y)] daiX) 



lim 



/[Ag.Ai]^ *^ [/J' V5(A, x) dz^(x) y) dr(y)] da{X) 
lA'fiK^y)-f{y)dy 



The assertion of the theorem follows immediately from (j30p once it is shown that 



(31) 

Observe that 



lim 



\iXt £ A) 



lim 



u[E-{[X^,Xi])e-'^''lA\ 



t^oo F^Xt < z) t^oo V [E^{[X^, Ai])e-*^'^1 



[0,^1 



V e 



-tL" 



u [i?«([0, AiDe-'^'^l^] + u [i?-((Ai,oo))e-*^''l^] 



(32) 



i/[E-([0,Ai])e-*^n^] 



1 + 



z.[i?«([0,Ai])e-*^n^] 
i/[E-((Ai,oo))e- 



[E-{%Xi])e~tL^lA] • 



Since e *^''1a and £'^([0, Ai])e *^"1a are continuous, the function £"*((Ai, oo))e must 



-tL 



tL" 



also be continuous. Thus z^[-E'^((Ai, oo))e 1a] is well defined. By Lemma l2.1 



(33) 



1 



lim -\og\{u,E^{{Xi,<x,))e-''^^lA} 



> Ai. 
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As (/?(Ag,x) > for every x £ (0,oo) there is, by continuity, Ai > Aq such that for every 
AG [Ag,Ai] 

1 (p{\,x)dv{x) and / V9(A, y)7(7/)(iy are both positive. 
Jo J A 

Then it is easy to see that 

(34) - hm - log [ e'^' H ^{\x)du{x) [ ^{X,y) jiy)dyda{X) < Ag. 

t-^oo t J[A5,Ai] Jo J A 

Equation ([32]) . ([33]) and ([3l|) combine to show 



lim - J. 

t^oo jy[£;«^([0,Ai])e-*-^"lA] 

and therefore (ISTI). □ 



3.6. Entrance boundary at oo. As mentioned above, with exception of the recent [9], 
work on these problems has generahy assumed that is regular and oo natural. Intuitively, 
we should expect the problems to be easier if oo is an entrance boundary. We show that this 
is indeed the case in Theorem 13.161 as the spectrum of the operator L** is purely discrete. 
This interesting fact has not been mentioned by previous authors (c/. section 3 in [9]) working 
on quasistationary distributions for one-dimensional diffusions. The proof relies on standard 
ideas from the spectral theory of differential operators. 

Theorem 3.16. //oo is an entrance boundary, then the spectrum of L'^ is discrete. 

Proof. Assume that is a regular boundary point, and we begin by considering the case 
K = 0. Let / be a solution to the eigenvalue equation —^(7/')' = Xf on (0, 00), for some 
A > 0. 

Let x > 1 be any local maximum, and x > x the first local minimum following x (assuming 
there is one). Since / solves the equation tu = Xu one easily sees that local maxima are 
positive and local minima negative. Integrating by parts and using the fact that f'{x) = 
f'{x) = 0, we have 



0</(x)-/(x)= [ ny)dy 

J X 

\i{y)f'{y))i{y)-^dy 



\i{y)ny))' j\{z)-^dzdy 
-2A [%{y)f{y) f%{z)-'dzdy 



<2X{f{x)-f{x)) / 7(y) / lizY^dzdy 



from which we conclude that 



^ < / 7(y) / ^dzdy. 
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Since we have assumed that oo is an entrance boundary, we know that 



oo> 7(y) / 7(2;) dzdy 



, -s J X Jl 

pairs (x,x) 

> (2A)~^ • #pairs 




Since the zeroes of / are separated by alternating local minima and maxima, it follows that 
/ has only finitely many zeroes on (l,oo), hence only finitely many zeroes in all. It follows 
from a theorem of P. Hartmann [52\ Theorem 1.1] that the spectrum of (the operator with 
K = 0) is discrete. 

Suppose now that the spectrum of L'^ is not discrete. Then there is a such that Ex, has 
infinite-dimensional range. All such x are in the domain of and satisfy q'^{x,x) < A*||x|p. 
But then they are also in the domain of and satisfy q^{x,x) < A*||x|p. By the minimax 
principle for the discrete spectrum (c/. |53i Theorem 8.8]), this contradicts the fact that L 
has discrete spectrum. □ 

Remark 3.17. There are general necessary and sufficient conditions for the discreteness of 
the spectrum of Sturm- Liouville operators obtained in [13j, of which Theorem \ 3.16\ may be 
seen as a special case. However, general versions found in the literature, such as the main 
result in [13] and Theorem 1 in [41j, do not seem to be immediately applicable. 



In this section we consider the problem of convergence to the Yaglom limit. More precisely 
we ask for conditions, which ensure that Xt converges to the quasistationary distribution 
given by the density '^{Xq, •). Recall that we always assume that is regular. Except in the 
final subsection 14.51 we will also assume that 00 is a natural boundary point. 

4.1. High killing at 00. In this section we consider the case where the asymptotic killing rate 
is strictly bigger than Aq. Theorem l2.6l shows that one has convergence to the quasistationary 
distribution if and only if the lowest eigenfunction is integrable. We give a new proof of this 
assertion and moreover prove that the lowest eigenfunction is actually always integrable. 
Therefore lim inf k > Aq always implies convergence to the quasistationary distribution. In 
contrast to j47j . we do not need to assume that 00 is a natural boundary: Since the bottom of 
the spectrum is an isolated eigenvalue, the corresponding eigenfunction is square-integrable, 
as well as Ag -invariant. 

Theorem 4.1. Suppose that liminfa;_!.oo '^(a^) > Aq. Then we have 



where uxk. G 2? denotes the uniquely determined (up to positive multiples) eigenfunction 
associated to the eigenvalue Aq. Moreover, the process {Xt) associated to the Dirichlet form 
q'^ converges to the quasistationary distribution . 

The theorem will be the direct consequence of two lemmas: Lemma 14.2^ which states that 
quasilimiting convergence follows whenever the eigenfunction ti^g is in S? and 2} ; and Lemma 
14. 3|, which states that is indeed in £} when lim'vaix^oo i^ix) > Aq. 



4. Convergence to Quasistationarity 



(35) 
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Lemma 4.2. Suppose G -G^ n Then (l35|) holds. 

Proof. We know from Lemma 13.31 (part |vj) that Ag is an isolated eigenvalue. Therefore, the 
eigenfunction ux^ is square integrable and satisfies 

(36) e~*^ = e~^^oux^ in hence identically (since ux^ is continuous). 

By (I22p . for r > sufficiently small, 



P''(i,a;,y) 



<C{x 
<Cix 



IbAx) uxs{x)l{x)dx 
IbAx) P'^i* + 1' ^' {x) 7(x)da 
/b.(x)^A5(x)7(^)c?^ 
e-(*+i)\"uA5(y) 



ux^{x)-/{i)dx 

For fixed x, p'^{t,x,y)e''^o is dominated by a constant times ux^{y), which is in £,^. The 
dominated convergence theorem, together with (120p . implies that there is a constant c such 
that for any Borel set U, 



oo 

e^5^p''{t,x,y)lu{y)7iy)dy 



r 

lim e^o*Px.(Xt GU,Ta>t) = lim / 

= cuAg(a;) / ux^{y)-i{y)dy. 



Taking quotients, 



lim P^(Xt G [/ ra > t) = lim — ^— ^ '- 

_ limf_,oo eW^(Xf e^,Ta>t) 
limt^ooe^o*P^(r« > t) 

_ c/f;UAg(y)7(y)c^y 



□ 



For the second part of the proof we apply an argument used in [8] to derive properties 
of eigenfunctions of Schrodinger operators. Some modification is required to deal with the 
complication that we have a domain with boundary, and we do not know a priori that 
the eigenfunctions are bounded. The one-dimensional setting helps us to overcome these 
complications. 

Lemma 4.3. Assume that Xq < K := liminf^^-^oo ^(2;)- Then the square integrable non- 
negative eigenfunction ux^ associated to the isolated eigenvalue Aq is integrable with respect 
to the measure T. 



Proof. By (I36p and the Feynman-Kac formula 



(38) e-^oiux,{x)=E, 



■^o-(^s)ds^^^^(^Xt),To 



> t 
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for every x £ [0, oo), where X* is the diffusion which is killed only at the boundary. For t > 
we define the martingale 

By the assumption Xq < K there exist positive real numbers a and e such that k{x) — Xq > e 
for every x £ [a,oo). Let Ta be the first hitting time of the set [0,a]. 
By the optional sampling theorem we get for every T > and x > a 



^Ag (^TaAT)l{To>T,AT} 



(39) 



< e-^'^E4nAg(X^)l{To>T}] +nA5(a)E, 



-eTaAT 



By Lemma 12.11 and the Spectral Theorem ([5]) the first term is bounded by 



oo \ 1/2 

2TA rii pO^.. ||2/ 



C^{x){ 1^ Xe-''^d\\E%xsr{X) 



(40) 



oo \ 1/2' 

2TA.||TpO„. ||2 



<e-^^2T-i/2| 



T-^oo 



0. 



(41) 



We have then, from (j39p and the Dominated Convergence Theorem, that 

eTaATj 

UA5(a)E. [e-'^'^] 



< -u-AsCa;) < lim nA-(a)Ea; le' 

" T— i-oo " 



We now appeal to a basic fact from potential theory (stated and proved in much greater 
generality as Proposition D.15 of [T6]; see also p. 285 of [6]): There is a constant C{a) such 
that for all X > a, 

(42) E4e-^'^»;r, < oo] = C{a)g'{x,a), 

where g"^ is the e-potential, defined by 

("OO 

/(a;,y) = / e-^*p{t,x,y) dt, 
Jo 
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where p{t,x,y) = p^{t,x,y) denotes the integral kernel of the operator e Since e *^ is 
self-adjoint, the integral kernel p{t, x, y) is symmetric with respect to F, so that 



J a 

POO 

<C{a)ux^^{a) g^{x,a)'y{x)dx 





fOO f'OO 

C(a)uAg(o) / / e~'^^p{t,x,a) 'j{x)dx dt 



C{a)ux^ (a) / / e '^^p{t,a,x) j{x)dx dt 
Jo Jo 



< C{a)e ^ux^{a). 

Since ux^{x)'y{x) is bounded on [0, a], this completes the proof. □ 

Remark 4.4. The above result reflects a general principle, which seems to be well-known 
to analysts and mathematical physicists: The decay of the eigenf unctions associated with 
isolated eigenvalues is dictated by the decay of the Green's function, at least in regions where 
the potential k is negligible . 

4.2. Low killing at oo: The recurrent case. We assume for the remainder of this section 
that K := lim^^^oo k{x) exists. Whereas the total surviving mass in the case -ftT > Aq declines 
at the strictly exponential rate e~^o*, in the case limjj^oo i^{x) < Aq one typically has 

(43) lim e^o* F^{Xt G ^, > = 

t— >oo 

for every bounded Borel set A C [0,oo). (This can be seen for a Brownian motion with 
constant drift by direct computation.) Equation (I43p remains true for every diffusion, if the 
bottom of the spectrum of the diffusion generator is not an eigenvalue in the £^-sense. Thus 
we cannot rely upon arguments that assume a spectral gap. 

It may seem surprising that, despite the complicated relationship between the unkilled mo- 
tion and killing for determining the lifetime of the process (and hence, whether it returns to 
its starting point), the conventional transience/recurrence dichotomy for the unkilled process 
is exactly the criterion that distinguishes between convergence and escape to infinity. We 
begin in this section by assuming that the unkilled process is recurrent, which is equivalent 
to assuming that '^{x)^'^ dx = oo, and show that this implies convergence to quasista- 
tionarity. In particular the lowest eigenfunction 99(Aq,-) is integrable (but now not square 
integrable) with respect to T. In section 14.31 we then address the case when the unkilled 
process is transient. 

As in |47] we define, for Borel sets A, the family of measures 

Ft{v,A)=¥,[Xte A\TQ>t) 

and 

at{v,r) = F^{Td>t + r\ Ta>t)= j Ft{i^,dy)¥y{TQ > r). 

If the process Xt started from the compactly supported initial distribution escapes to in- 
finity, then for any sequence (in)nGN converging to infinity the measures Fj„(z/, •) converge 
weakly to point the niGcisurG Sqq. If the process Xi started, from u converges to the qiia>si- 
stationary disdistribution ip then then the hmit of Ft^{u, dy) is concentrated on M+, and has 
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the density -i^ — y (^o \\ with respect to F. The next lemma is in essence a combination of 

Jo f\.^o'y)'y'yy)'^y 

Lemma 5.3 and Theorem 3.3 in |47], together with our Lemma 13.71 

Lemma 4.5. Assume that oo is a natural boundary point and suppose that Aq ^ K. Then 
the limit a{i^,r) = limt^ooCLtii^ir) exists, and satisfies 



(44) aiu,r) = F{iy,R+) ip{X^,y)¥y{re > r) -f{y)dy + {l - F{iy,R+)) 



Either F{i',M^) = for every compactly supported initial distribution v or = 1 for 

every such v. 

There exists rjy ^M. (called the asymptotic mortality rate) such that 
(45) a{v,r) = e-'^^\ 

If the process escapes to infinity then rji, = K . 

Proof. Let v he a, compactly supported initial distribution. Let (t„) C (0, oo) be a sequence 
converging to infinity. On the compactification [0, oo] of (0, oo) the sequence of measures 
Ft^[i',dy) has a limit point. By Theorem 12.61 this limit point is either a measure on (0, oo) 
which has the density -r^ — y(^o \\ with respect to the measure T or is the point mass at oo. 

Theorem 12.61 shows that there is only one limit point and that the limit point is independent 
of the sequence (t„) and the initial distribution u. Thus Ftiv^dy) converges weakly. If oo is 
natural, then 

lim ¥Jtq >r)= e-^''. 
y^oo ^ ' 

This shows that 

^hm I Ft{u,dy)¥y{Te > r) = F{u,m+) J ^{X^,y)Fy{Ta > r) jiy)dy + (l - F(i/, M+)e-^0' 



which is (|44|) . 

For any r, s > we have 



. , . r FA^d>t + r + s 
a(z/, r + s) = nm 



t^oc P(ro > tn) 

,. K{TQ>t + r + s)Fy{Ta>t + s) 
= hm ; — 

t^oc Fy(TO>t + s) F^{Td>t) 

^. Fy{T9>n + r + s) F^{T9>n + s) 
= hm — lim — 

n-^ca Fy{T9 > n + s) n^oo Fy[TQ > n) 

= a{v, r)a{v, s), 

which directly implies ()45p . The final statement follows directly from (j44p . □ 

In order to decide, whether Xt converges to the quasistationary distribution, we investigate 
the asymptotic behavior of the function r i— )• Fy{TQ > r) , as r — )• oo. 

Lemma 4.6. Suppose that the lowest eigenvalue Aq is strictly positive and that T is a finite 
measure. Then for any compactly supported initial distribution v we have 

- lim \ \ogF,{TQ>t) =Ag 
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Proof. Since F is finite, the constant function 1 is in Lemma |2 . 1 1 implies 



(46) 



1 1 /" K 

lim sup - log (rg > t) = lim sup - log / e^*'^ l{x)di'{x) 



<-K- 



Now we need a corresponding lower bound. Fix any z > 0, and let I := {x : \x — z\ < ^ A|}. 
By (p2]) there is a constant C(z) = sup^g2C('"^) such that 



|e-t^^lir = (e-*^''lx,lx) 



(47) 



p^it,x,y)driy)dr{x) 



IX JX 

<Ciz) l^p^{t + l,z,y)j{y)dy 
<Ciz)F,{T9>t + l). 
By Lemma 13.141 and Lemma 12.21 we see that 

liminf ^logP^ra > t) > - inf suppd||£^''lxf (A) = -A^, 

Since the exponential rate of decay of Fz{Td > i) is locally constant in z, by Proposition 13.13] 
this completes the proof. □ 

Theorem 4.7. Suppose that K < Ag, and j{x)^^ dx = oo. Then Xt started from an 
arbitrary compactly supported initial distribution v converges to the quasistationary distribu- 
tion with T -density proportional to (/?(Aq,-). Moreover, the asymptotic mortality rate rjy is 
independent of v and equals Ag . 

Proof. If Xt escapes to infinity then we know from Lemma 14.51 that 

F^jTQ >t + r) _Kr 

a[u,r) = lim — — r — = e . 

Since by assumption Ag > K, when a > part (|viip of Lemma 13.31 tells us that Aq > 0. The 
strict positivity of Aq together with the assumption 'y{x)~^ dx = oo allows us to apply 
part dni) of Lemma 13. 3| to conclude that the speed measure F is finite. When a = and F 
is infinite the same reasoning holds, leading to a contradiction. Therefore we may assume, in 
any case, that F is finite 

Therefore Lemma 14.61 shows that for every compactly supported measure ly 

- lim ilogP,(Ta>t) =AS. 

t— J-oo t 



In the case of escape to infinity equation (j45jl implies 

- lim I logP^(ra >t) = r]^ = K^X^. 

Therefore the assumption F(iy,M^) = cannot be true and thus by Theorem 12.61 we conclude 
F(i^, M+) = 1 and -F(z^, oo) = 0. Thus Xt converges from every compactly supported initial 
distribution to the quasistationary distribution V'(Aq, •). □ 
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The above theorem has the following Corollary, which in a slightly more restrictive form 
already appears in the work [12] of Collet, Martinez and San Martin. The proof presented in 
|12] suffers from a gap, so it seems to be worth presenting an alternative (and more general) 
proof of the assertion. 

Corollary 4.8. Suppose k = and oo is a natural boundary point, and the process Xt is 
recurrent, with a > 0. Then 

• if Xq > then Xt converges from every compactly supported initial distribution v to 
the quasistationary distribution with T-density proportional to (^(Aq,-)/ 

• if \q = then Xt started from v escapes to infinity. 

Proof. The first part of the assertion follows directly from Theorem 14.71 In order to prove 
the second assertion, observe that the function 

1 a fy , 
R{y) ■■= -r— + / 

1 + a 1 + a Jo 

satisfies LR = 0; since i?(0) = 1/(1 + a) and i?'(0) = a/{l + a) it is in Pq, and is the unique 
eigenfunction ip{0,y). Since Aq = 0, we have 

^{>^o,y)l{y)dy = — — / l{y)dy + — — / l{y) j 7 {x)dxdy = 00, 
1 + a Jo 1 + a Jo Jo 

by the assumption that 00 is a natural boundary. □ 
4.3. Low killing at infinity: The transient case. 

Theorem 4.9. Suppose that co is a natural boundary point and that j{x)~^ dx < 00. 
If K < Aq then Xt escapes to infinity from every initial distribution. The rate of escape is 
exponential with rate AfJ — K, in the sense that for all z,x > 0, 

(48) lim sup ^logr^Xt < z \ T9 > t) = -{X^ - K). 

Proof. Observe that the condition j{x)^^ dx < 00 implies that for each a G (0, 00) and 
each X G (a, 00) the unkilled diffusion (corresponding to the generator L) started from x has 
nonzero probability of never hitting a. For e > we can choose a = G (0, cxd) such that 
k{x) G {K — e,K + e) for every x G [a, 00). Then we have for every x G (a, c«) 

Px(t9 >t)=E, e-^'"(^^)'^^ To > t 

>e-(^+^)¥4r, = oo) 

Since Px(Ja = 00) is an increasing function of x, we can apply the Markov property to see 
that there is a nonzero increasing function C{x) such that for all x > 

(50) P^(Ta >t) >F^{Xi>a + l) ■ ^mf^^P^./(Te > t - l) > C7(x)e~(^+")* 

On the other hand, for any fixed z > we can apply the bound (I13p and Lemma 12.21 to see 
that 

M^t <z,Ta>t) = (e-*^''l[o,,])(x) 

<(C„(x)AS + C;)||l[o,.]||e-*^S 
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for all t > 1/2 Aq. Combining (|50p and (|5ip. we see that there is a constant C such that 

c 

We conclude that for all x > 0=-, 



(52) M^t <z\T9>t)< ||lfo,,]||^e-(^5-^-^)*. 



limsup - logF^{Xt <z\Ta>t)< -(AfJ - K) + e. 

By Proposition 13 . 1 3l since e is arbitrary, we conclude that the limsup is no more than — (Aq — 
K). 

In particular, we have shown that the process escapes to infinity. By Lemma [4. 51 it follows 
that limj_^oo ^x{Td > t + 1 \ Tg > = , from which follows that 

lim t-'^logFJrs > t] = -K. 

Lemma 13.141 tells us that 

lim r^\og¥.,{Xt<z} = -X^^, 

from which (j48p follows immediately. □ 

If K is eventually constant — that is, for some a we have k[x) = K for all x > a — then 
we can strengthen the conclusion of Theorem 14.91 slightly. 

Corollary 4.10. Suppose that k is eventually constant and that Aq > 0. Then for every 
x,z £ (0, 00) 

sup 

Proof. If K is eventually constant, then (j49p and ()50p hold with e = 0, hence (|52p as well. □ 

Remark 4.11. The case k = corresponds to the setting considered in [35]. Theorem 4 of 
\35\ includes a slightly weaker version of the result in Corollary \4-10\ obtained by different 
methods. The above theorem shows that when k = the S!^ -eigenvalue Aq gives the exponential 
convergence rate at which Xt escapes to infinity. 

As already mentioned in Remark 12.51 a quasilimiting distribution i), which in our case is a 
probability measure on (0, 00) is always quasistationary in the sense that for every Borel set 
A C (0,00) 

Pi>(Xt £A\T9>t)= v{A)- 

but the converse need not hold true. In the cases where we know there is no quasilimiting 
distribution, though, because the process escapes to cxd, we can show that there is also no 
quasistationary distribution. 

Corollary 4.12. Let 00 be a natural boundary, with Xq > K and ^{x)~^ dx < 00. Then 
there is no quasistationary distribution. 

Proof. Assume that P is a general quasistationary distribution. The measure D is absolutely 
continuous with respect to F with a positive continuous density g : [0, 00) — )• (0, 00) (for a 
sketch of the proof of this fact we refer to the appendix). There is a A such that Pj>(Tg > t) = 
e~^* . By (f50]l . for any positive e, 

(53) e-^* = P^jra > t] > e-(-^+^)* j C{x)du{x), 
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which means that X < K. For any fixed xq > 0, 
Pr>{Xt<^, ra>t} = (5,e-*^^l(o,,]> 

(54) = (5l[o,a;o]'e~*^ 1(0,^]) + (^l{xo,oo),e"*^ 1(0,^]) 

< ||5l[o,a;o]|| • lle^'^^lcc^lll + sup ¥^{Xt <Z,Te> t} 

X>Xq 

Since ||(7l[o,a:o] II ^^^d ||l(o,^]|| are both finite, we can use ([5]) and (|5T]l to see that there is a 
constant B such that 

(55) P£>{Xt < z, ra > t) < 5e-*^o. 
Combining ([53]) and (f55|) . we see that for all positive t, 

(56) i>([0, z]) = ^,{Xt <z\Te>t]< Be-'-^o-K)t^ 

so D must be identically on [0,oo). □ 

4.4. Processes that may not hit 0. Consider a process which is killed only at (that 
is, with K = 0.) If the process is not almost surely absorbed at eventually — that is, if 
Pa; (To = oo) > — we may wish to condition the process at time t on being killed eventually, 
but not yet. That is, we consider the long-time asymptotics of 

^x{Xt G • I To G (t,oo)). 

Conditions of this kind can often be found in the analogous problems in the theory of branching 
processes. This problem can be reduced to our previous analysis by an h-transform. The 
function h{x) = W^iTo < oo) is harmonic, and by general theory (see |40) . chapter 4, sections 
3 and 10) the process (Xt) conditioned to hit corresponds to the generator whose action 
is given by 

lV = (iT(/./))(x) = -^f"{x) + (-6(x) - ^)/'(^) 

The process associated to the operator can again be defined by Dirichlet form techniques, 
and the associated family of measures on the path space is denoted by P^; . As explained above 
we have 

P,(- I To < oo) = P,(-). 

The operator can be realized as a self-adjoint operator on the Hilbert space -C^((0, oo), h{x)'^'y{x)dx) 
The transformation V : £'^{{0,oo),h{x)^^{x)dx) — ^ £^((0, oo), 7(x)(ix) defined by Vf = fh 
is unitary, and defines a unitary equivalence between L and L^, so the spectrum is invariant 
under /i-transforms. In particular, positivity of the bottom of the spectrum of L implies the 
positivity of the spectrum of L^. Since absorption is certain with respect to the measure P^; 
we can apply our previous results in order to conclude that for every Borel set Ac (0, oo) 

, (p^{\n,x)h{X)^(x)dx 

lim ¥^{Xt e A I To e (t,oo)) = ;°oo:, °' ')J,) ' , 
t^oo Jq ip"'[Xo,x)h[X)^[x)ax 

where ip^{XQ,x) is the unique solution of {L^ — Xo)u = 0, which satisfies (p^{X(j,0) = and 
(<^'^)'(Ao,0) = l. 
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Remark 4.13. It seems to be a rather general principle that there are three possibilities. 
The first possibility is the non-existence of quasistationary distributions. If there exists a 
quasistationary distribution then it is either unique or there is a whole continuum of quasi- 
stationary distributions parameterized by a real interval. This is at least true for birth and 
death processes on the non-negative integers, cf. [IQj. 

4.5. The Case of an Entrance Boundary at oo. Observe that 'j{x)~^dx = oo if oo 

is an entrance boundary. This follows from the fact that in this situation the total speed 
measure Jq ^{x)dx must be finite. Thus, the situation is essentially the same as in Theorem 
14.71 Indeed, we always have convergence to quasistationarity if oo is an entrance boundary. 

Theorem 4.14. Assume that is regular and that oo is an entrance boundary. Then the 
bottom of the spectrum is an isolated eigenvalue with associated non-negative eigenfunction 
u\i^. From every compactly supported initial distribution v, the process Xf converges to the 
distribution with density ux^ / u\^{x)'y{x)dx with respect to T. 

Proof. The first assertion follows from Theorem 13.161 Lemma 14.21 directly implies that Xt 
converges to the quasistationary distribution u^g from every compactly supported initial 
distribution if and only if ux^iu) l{v)dy is finite. Since we are assuming that is regular 
and oo is an entrance boundary the speed measure F must be finite. Thus the 2?{T) function 
ux-^ is also in £-'^(r). □ 

Remark 4.15. The uniqueness of quasistationary distributions (assuming their existence) 
was addressed in the recent paper [9J, for the case k = 0. There it is shown that the uniqueness 
of quasistationary distributions in this case is equivalent to the assertion that for any a > 
there exists ya > such that sup^^^^ E^; [e""^'''' ] < oo, where Ty^ denotes the first hitting time 
of ya- Thus uniqueness of quasistationary distributions is equivalent to the 'time of implosion 
from infinity into the interior' having exponential moments of all orders. It is also proved 
in [9j that this is equivalent to oo being an entrance boundary. If k = 0, and if absorption 
is certain, our results show that the existence of a quasistationary distribution is equivalent 
to the existence of exponential moments of the first hitting time of 0. Both results together 
account for the existence and uniqueness of quasistationary distributions. 

Appendix 

In this Appendix we sketch a proof of the regularity of quasistationary distributions of 
one-dimensional diffusions with one regular boundary. 

Lemma A. Let L'^ be one of the selfadjoint realizations considered in this work, of the 
Sturm- Liouville expression r + k in SP; and let D be a quasistationary distribution. Then v 
is absolutely continuous with respect to the measure T, with a positive and continuous density 
g : [0,oo) ^ M. 

Proof. The main assertion of the Lemma will be almost obvious to readers who are familiar 
with regularity theory for stationary distributions. Indeed, the main strategy we follow is 
very similar to the case of stationary distributions. Straightforward arguments show that D 
is absolutely continuous with respect to the measure F. Denote by g the density of D with 
respect to F. The equation 

e-^*z>(/) = E^[f{Xty,Ta > i]; A > 0, / G a((0,oo)). 
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which results from quasistationarity of u, imphes that 

(57) V/ G Cr ((0, oo)) : {D, (L'^ + A)/) = J g{x){L^ + A)/(x) dr(x) = 0. 

This means that for any < c < d, 

g G P(r*J and T^^g = 0, 

where T*^ denotes the adjoint (taken in the Hilbert space £^((c, d), F)) of the minimal operator 
Tc^d defined as the restriction of the differential operator L'' + A to C^((0, oo)). The domain 
of T*j^ is given by 



'D{T*^) = {/ G £^((c, (i),r) \f,"ff' absolutely continuous in {c,d) and 

{lf')' + {K + \)f€S.\{c,d),T)} 



(58) -1, , . , 
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and for / G P(T*^) one has T*^ = 27(7/')' + (/^ + A)/. Since c, d G (0, 00) are arbitrary we 
conclude that 

(59) ^(7/)' (^) + (^ + A)/ (x) = in (0, 00) 

Due to the regularity of the boundary point we conclude (using standard ODE theory) that 

lim r(x) G M 

exists. □ 
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